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A real-time functional-integral method is used to derive an effective ac- 
tion that gives the transmission spectrum of a tunneling particle interacting 
with a bath of harmonic oscillators. The transmission spectum is expressed in 
terms of double functional integrals with respect to the coordinate of the par- 
0> ' tide which are evaluated by means of stationary-phase approximation. The 

equations of motion for the stationary-phase trajectories are solved exactly for 
an arbitrary spectral density function of the bath, and the obtained solutions 
are used to find the transmission spectra for specific examples. For a bath 
I with single frequency u>, an analytic expression of the transmission spectrum 

is obtained which covers from sudden tunneling {ujTq ^ 1) to adiabatic one 
^ {loTq ^ 1), where Tq is the time it would take a classical particle to traverse 

Q>^ ■ the inverted bare potential barrier. For a bath with Ohmic spectrum, the 

. differential tunneling conductance at low bias voltage V and for ryTg <C 1 

is found to obey a power law ~ {eVTo/n)'^'^°^o/'^''^, where r] is the friction 
' coefficient and 5*0 is the tunneling exponent in the absence of interaction. 

0\ '. PACS numbers: 73.40.Gk 



B 

o 
o 



I. INTRODUCTION 



Persistent efforts have been devoted to the development of computational methods in 
tunneling. The instanton method offers a way of evaluating the decay rate of a metastable 
state |ll],|2|. The central idea of this method is to deduce the information on the lifetime 
of a metastable state from the imaginary part of its free energy that could only be defined 
! by means of analytic continuation. The instanton method, however, does not give the 
transmission spectrum — the energy distribution of transmitted particles — because the free 
energy is calculated from the partition function which is defined without reference to the 
final states of the transmitted particles. A modified version of the instanton method which 
takes only half of the closed-loop bounce trajectory can compute the tunneling rate for the 
case in which both initial and final states of the bath are in the ground state |Q, but the 
instanton method has not succeeded in finding full transmission spectrum ||^. 

Quite often, however, one encounters situations in which one needs to find the transmis- 
sion spectrum. Examples are the Coulomb blockade of tunneling deep inelastic collisions 
of heavy ions 0, atom-surface scattering [0, chemical reactions [H, and, of course, tunnel- 
ing spectroscopy Ueda and Ando have recently computed full transmission spectrum 
of a tunneling particle interacting with a bath of harmonic oscillators by combining an 
operator-algebraic method and a functional-integral one followed by analytic continuation 



to imaginary time [10]. Real-time approaches have also been used to find information on 
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transmission spectrum. Mohring and Smilansky used the influence functional to 
calculate the variance of the energy that a tunneling particle exchanges with the bath. By 
applying a time-dependent WKB method to approximate functional integrals for the coor- 
dinate of the particle and then by performing functional integrals for the bath degrees of 
freedom, Bruinsma and Bak derived the effective action to the lowest order in the coupling 
constant 0]. 

The primary purpose of this paper is to use a real-time functinal-integral method to 
derive an exact expression of the transmission spectrum of a tunneling particle interacting 
with a bath of harmonic oscillators, and to show how to get information on the tunnel- 
ing exponent from the exact expression by means of stationary-phase approximation. The 
method presented in this paper allows a systematic expansion of the effective action in terms 
of the coupling constant and takes account of the fact that the time T it would take the 
tunneling particle to traverse the barrier along the stationary-phase trajectory is modifed 
by the interaction — a point that has eluded the treatment using a time- dependent WKB 
approximation 0. The present method also allows us to analytically study the dynamics 
of tunneling from the sudden limit ujTq ^ 1 to the adiabatic one ujTq ^ 1, where uj is the 
characteristic frequency of the bath and Tq is the time it would take a classical particle to 
traverse the inverted bare potential barrier. 

This paper is organized as follows. Section |I| obtains a general formula for the trans- 
mission spectrum (Sec. [11 A]) and applies it to the Caldeira-Leggett model (Sec. [11 B|) ||T2 



where the transmission spectrum is expressed in terms of double functional integrals with 
respect to the coordinate of the particle alone. Section |T| evaluates the double functional 
integrals by means of stationary-phase approximation and exactly solves the equations of 
motion for the stationary-phase trajectories for an arbitrary spectral density function of the 
bath. The obtained solutions are used to analyze the case of a single-frequency bath in 
Sec. 13 and the case of an Ohmic bath in Sec. 0. The relation of the present method to 



the instanton method is discussed in Sec. [Ill D| . Section [IV (J[ examines the distribution of 
sideband intensities of particles transmitting through a quantum-mechanically fluctuating 
barrier and compares the results with those obtained for the case in which the barrier is 



time-modulated in a prescribed manner |T3| 



II. FORMULA FOR THE TRANSMISSION SPECTRUM 

A. General theory 

We consider a situation in which a particle of mass M tunnels through a potential 
barrier V{X) while interacting with other degrees of freedom - "the bath." A general 
scattering-theory formula for the rate of transition from state |\E'i) to state |\l'f) is given 
in the Schrodinger representation by 

Ti^f = lim -^|(vI/f|Texp(-^ f' Hdt) \^-:)\\ (1) 

where H is the total Hamiltonian and T is the time-ordering operator. When H does not 
depend on time explicitly, the factor Texpf— jl^Hdt) may be replaced by e-'^^i^^t^^O/^^ 
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When we are not interested in the state of the bath, we make a statistical average of Eq. (|l]) 
over initial states of the bath and sum over final states: 

r{E,E') = ^hm^^X:P(ei,)|(vl/f|Texp(-^£i7dt) m'6[E - E' - {el - e^)], (2) 

where E is the initial-state energy of the particle, ej-, and are the initial-state and final- 
state energies of the bath, and -P(ej^) denotes the probability distribution of ej_,. The delta 
function restricts the energy gain {if E — E' > or loss otherwise) of the bath between the 
initial state and the final state to be equal to E — E' which at the limit of tf — ti = oo becomes 
the energy loss ( if — £" > or gain otherwise) of the tunneling particle. Therefore E' 
has the meaning of the final-state energy of the particle and Eq. may be written as 

T{E, E') = ^^hm^ ^ X: I exp (-^ \E, e\,)\'5[E - E' ~ {ei - . (3) 

Inserting into Eq. the completeness relations for the coordinate of the particle in the 
initial and final states, we obtain 

T(E,E')= lim r ^e^^'^-'^'^'IdxJdxJdYJdY, 



where 



tf4->±oo — tiJ-oo27rh 

X {E'\Xf){Xi\E){E\Y{) {Yi\E')AiXi, X,, Yf, Y,; t), (4) 



A(Xf, Xi, Ff, Fi; t) = Yl P{el){Xi, e[\Texp f^^Hdi^ \X,, ej,) 

x{Y„e\,\fexp(^^£Hdt^ \Y,,ei) e"^(<-j>)*, (5) 



where T is the counter-time-ordering operator. Equation together with Eq. (^) gives 
the transmission spectrum that a particle with initial-state energy E transmits through a 
barrier with final-state energy E'. This formula applies whether or not the situation of our 
concern is related to tunneling. The condition that the particle indeed tunnels through the 
barrier will be imposed later [Eq. (pS])]. 



B. Effective action for the Caldeira-Leggett model 



By explicitly evaluating Eq. for the Caldeira-Leggett model ||T2[, we derive an ex- 
act effective action for the transmission spectrum of a particle interacting with a bath of 
harmonic oscillators. The Hamiltonian of the Caldeira-Leggett model is given by 

where X and P are the coordinate and momentum of the particle; Xa and Pa are the 
canonically-conjugate generalized coordinate and momentum of mode a of the bath, and 
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rric, and uja are its mass and frequency; and fa{X) is a function that characterizes the 
interaction between the particle and the bath. 

The Caldeira-Leggett model includes the counter term X]«("^ai^a/2)/a(^)- Whether or 
not this term should be included depends on the problem concerned, but we proceed with 
it as it can easily be subtracted at any stage in the following discussions. Because the 
transmission spectrum of a particle is meaningful only if the particle is free in the initial 
and final states, we require that 

V{X,^,)) = V{Y,^,)) = (7) 

and 

= = ffK (8) 

where f^^^ may depend on parameters such as rua, cOa, etc., but should not depend on the 
dynamical variables of the particle. If these requirements are met, the particle and the bath 
are decoupled in the initial and final states. The energy eigenstate of the bath in the initial 
(final) state \e^^^) is represented by with its eigenenergy given by 

ef = E hu^mf^ + l), (9) 



{mf>} 



where mj^^-* denotes the excitation number of mode a of the bath. 

To eliminate the bath degrees of freedom, we express Eq. (|]) in coordinate representation 
of the bath by inserting the completeness relations: 



/CO roc poo poc 

dXal / dXa2 / dXa3 / dXai 

X (Xf, {x^,}\e-^nm-tO\x,, {x^2}){Y,, {x„3}|e^^(*^-*')|Ff, {x^,}) 



X 



Here (x|m^) and {x\n^^) are given in terms of the wavefunction ipni^) of the harmonic 
oscillator as ipma(^^ ~ /a) ^-^id ipna(^^ ~ respectively. The summation over nl^ gives 
the Green's function of the shifted harmonic oscillator: 



^ y 27rih sin Uat \ 2ih sin Uat J ' 

where Xai{4) = 3^01(4) — fa- When the initial state of the bath is at thermal equilibrium with 
temperature (3~^, P{rn}^ is given by 

P(m'a) = (1 - e-^'*^-)e-'^'"-^'^", (12) 
and the summation over can be carried out to yield 
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2 sinh 



X exp 



-2171% siiiuja{t + i(3h) 



—2ih sin uja{t + i(3h) 



(13) 



where Xq2(3) = a;Q,2{3) — /„• The probabihty amphtudes (Xf, {xoijle r'^'-*'' {xa2}) and 

(Xi, {xQ,3}|e^^(*t~*i)|Xf, {xa4}) may be expressed in terms of functional integrals over the 
coordinate of the particle and over the coordinates of the bath. Because the Hamiltonian 
is quadratic with respect to the coordinates of the bath and the interaction term is linear 
in them, path integrations over the bath degrees of freedom can be carried out explicitly, 
giving PI 



(Xf,Ki}|e-i^(*'-*-)|Xi,K2}) = n^ 



2TTih sin Ua{t{ — ti) 



rX(t{)=X{ rt{ 

X // VXexpl- dti 

'lx{ti)=Xi [rUti 



^X\t^)-V[X{t,)] 2 



inrtnUOn 



2hsmuJa{ti - ti) 

^li cosuJa{tf-ti)-2xaiXa2 + 20;^^ dtifo, [X{ti)] sin uj^ti -t;) +Xa2 sin t^«(tf -ti) 

2uj^J^ dtil^dt2fa[X{ti)]fa[X{t2)]smUg{tf-ti)smUa{t2-ti) 



(14) 



and a similar expression can be found for {Yi, {xa3}\eT^"^^'~^'^\Yf, {xai}). Substituting 
Eqs. ([TT|), (|I^) and (|T^) into Eq. (|TD|) and performing integration over Xai {i = l,---,4) 
finally yield 



yl(Xf,Xi,rf,ri;t) 



X{ti)=Xi irY{tt)=Yt 

VX // VY exp 

X{ti)=Xi JlY[U)=Yi 



(15) 



where the effective action S'efr(t) is given in matrix form as 



5'eff(t) = dti 



fx^ - ^Y- - V{X) + V{Y) 



+ E 



— / dti I dt2 



4 sinh ^ 



X 



U[X{t,)],-^Y{t^)] 



coshwa^^ - « 1^1-^2!) cosha;„(^ + i{ti-t2-tj 
coshujai^ — i{ti—t2+tjj coshcUol^^ + i\ti — t2 



(16) 



This is the desired effective action for the Caldeira-Leggett model that gives the transmission 
spectrum of a tunneling particle interacting with a bath of harmonic oscillators. 



For the case of a separable interaction |]T2| in which fa{X) is given by fa{X) 
Cag{X) /maOjl^, Eq. (ITO) reduces to 



Ses{t)=£dh 



fx- - ft- - V{X) + V{Y) 



+t f'dh rdt2\g[X{ti)l -g[Y{ti)] 



a'{ti-t2) a<{ti-t2-t) 
a>{ti-t2 + t) a^"(ti-t2) 



^[^(^2) 



(17) 
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where the nonlocal integration kernels 



duJ{u)coshuj(^f-it 



2ti uj"^ 



sinh^ 



a>{t) = [a<{t)Y = 
a^{t) = [a^{t)Y = a>{\t\), 
are characterized by the spectral density function of the bath 



J(^) 



TT 



(18) 



(19) 



The Green's function matrix appearing in Eq. ( \17\ } reminds us of the closed-time-path 
Green's function (CTPGF) formalism |jl5[, which pursues the time evolution of the system 
by eliminating all pieces of information on the bath degrees of freedom. The present theory, 
by contrast, does not eliminate all of them but keeps the piece of information on how much 
energy has been transferred between the system and the bath. Hence we have the time 
argument t in the off-diagonal elements of the Green's function matrix in Eq. ([T^; it plays 
the role of giving the transmission spectrum after the Fourier transformation in Eq. (H). 

Substituting Eq. (^) into Eq. (^) gives 



TiE,E') 



lim — 

tf,i->±oo 



ti 



°° dt 



o2nh 



X 



fX{t()=Xf rrY{tf)=Yf 

VX VYexp 

X(ti)=Xi JJY{ti)=Yi 



e^(E-E')tJdX{JdX-,JdYiJdYi{E'\X{){X-,\E){E\Y^{Yi\E') 

(20) 



The integrals in Eq. (^) over Xf, Xj, Yf, give the boundary conditions on functional 
integrals 0. To see this, let us introduce the center-of-mass coordinate R= {X + Y)/2 and 
the relative coordinate q = X ~ Y . The relevant part of the effective action becomes 



dti 



^X'^—Y' 
2 2 



MRqdti = M R{t{)q{tf) - i?(ti)g(ti) - / MRqdti. (21) 



Since there is neither the potential nor the interaction in the initial and final states as 
required in Eqs. (|^ and (|]), there is no acceleration R. Noting that 



{E'\Xi) = e 



-ikfXf 



(22) 



where h = \l2ME/h^ and h = ^2ME'/h, we find that 

ldX,JdX,JdYfJdY,{E'\Xi){X,\E){E\Y;){Y,\E')e^n[M^^^^^^ 

= {2TihY [5{MRi - hh)S{MRi - hk)f . 
This provides the boundary conditions 



(23) 



MR{tf) = hkf = V2mE', MR{U) = hh = V2mE (24) 
on functional integrals in Eq. (120). With a proviso that these conditions are met, we obtain 
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T{E,E')= lim 



1 



dt 



eUE-E')t 



X(ti)=Xi i..Y(ti)=Yf 

VX VYexp 

X{ti)=Xi JJY{ti)=Yi 



(25) 



This result together with Eq. (|TBp or (|T^ gives the desired exact expression of the transmis- 
sion spectrum of a particle interacting with a bath of harmonic oscillators. 

The integral of r{E,E') over E' gives the total tunneling rate TtotBi{E) for a particle 
with initial-state energy E. From Eq. (pS]) we obtain 



total 



T{E,E')dE' 



lim 

if.i^ioo tf — ti 



Xiti)=Xi irY(tt)=Yt 

VX PFexp 

X{ti)=Xi MY{ti)=Yi 



h 



Ses{t = 0) 



(26) 



Thus the total tunneling rate is given by the value of Ses{t) at t = 0. This action should 
be related to the imaginary-time effective action derived by Caldeira and Leggett using the 
instanton method because the latter also gives the total tunneling rate for the same system. 



We will show this in Sec. Ill D 



The physics contained in the effective action 

5eff(t = 0) + [S'efr(t)-S'eff(t = 0)]: 



is revealed if we rewrite it as Sesit) 



l-tf l-tf l-C 

S,s{t) = S,s{t = 0)+i dti / dt2 

Jti Jti Jo 



\lg[X{1:r)]g[Y{t2) 

1-^ 



+nB(cu)(l-e*"*)e 



iujt\ —iuj(ti-t2) 



(27) 



where n-Q{uo) = (e^^'^ — 1)~^ is the boson occupation number. As shown in Eq. ( p6D the 
first term Sesit = 0) on the right-hand side of Eq. (|27|) describes the total tunneling rate. 
The second term comprises two ingredients in the curly brackets: the first of them describes 
emission of the energy quanta and the second describes absorption of energy quanta. At zero 
temperature, the second one vanishes because the bath has no energy to give off, whereas 
the first one survives because of the presence of spontaneous emission. The presence of 
spontaneous emission is a unique feature of the interaction with dynamical fields and makes 
a distinction from the interaction with classical fields as will be illustrated in Sec. IIV C . 



III. STATIONARY-PHASE APPROXIMATION 
A. Equations of motion 



We evaluate the functional integrals in ( pS]) by means of stationary-phase approximation 
in which the functional integrals are replaced by ordinary integrals along the stationary-phase 
trajectories. Assuming that the particle interacts with the bath only under the barrier, we 
choose 



UX) 



ifX<0; 
^X iiO<X <d: 



niauj^-- (28) 

-^d if X > (i. 



where and d are the left and right turning points of the barrier. Since the particle must be 
on one side (or the other side) of the barrier in the initial state (or the final state), fa{X) 
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given in Eq. (pS]) indeed satisfies the condition (||). The first condition (0) is also satisfied 
by our choice of the potential described later. 

The fact that the Hamiltonian does not depend explicitly on time means that the sys- 
tem is invariant under translation of time. The quantum fiuctuations around the stationary- 
phase trajectories therefore include the zero-frequency mode which produces the prefactor 
proportional to tf — ti Contributions of the remaining quantum fiuctuations produces 

an extra prefactor N which, in general, depends on energy, Tx, and Ty, etc. However, we 
assume that in the semiclassical limit its dependence is weaker than that of the tunneling 
exponent, and we will neglect its dependence in the following discussions. We thus find from 



Eq. @ that 



T{E,E') = N J 



dt 



.(E-E')t^l.S(€) 



(29) 



where S(t) is obtained from S'efr(t) by replacing the upper limits of integration by Tx (or Ty) 
for trajectory X (or Y) and the lower limits of integration by 0. The Tx and Ty will later be 
determined so as to meet conditions X{Tx) = Y{Ty) = d, provided that X{0) = Y{0) = 0, 
and may be regarded as (complex) traversal times in the sense that they are the "times" it 
takes the particle to traverse the barrier along the stationary-phase trajectories X and Y. 
At zero temperature S(t) becomes 



S{t) = 

Jo 



2" 



-X'-—Y^-V{X) + ViY)-^X^ + ^Y^ 



■ f 1 [ 1 f°°duj ^. 
+ 1 dti dt2 TT-^V^ 
Jo Jo Jo 2tt 



-'"i*i-*^ix(ti)x(t2)+e'"i*^-*^ir(ti)r(t2) -2e'"[(*^-^^)-(*^-^^)]x(ti)r(t2) 



+ id 



°duj J{u) ^ 
vr 



1-e 



'iuit\ 



Jo d 



Jo 



d 



(30) 



where it is understood that T is equal to Tx (or Ty) if the trajectory referred to is X 
Y), and 



.or 



E 



2 

TC Jo 



oo 



■doj. 



ijj 



(31) 



In Eq. ( pOD we have replaced the factor ^^^{ix-ii) in front of X{t])Y (t'2) by 
^iu)\{tx-Tx)~{ti-TY)\ _ This apparently ad hoc replacement must be made in order to get phys- 
ically reasonable results in the limit of cjTo 00 (see Sec. |^). In fact, in the limit of 
large loT^^ the only effect of the interaction between the system and the bath is the potential 
renormalization that should be canceled by the counter term. Thus there should be no effect 
of the interaction on the tunneling rate, which will be confirmed by the above replacement. 
If we could perform the functional integrals exactly, the original factor e*'^^*!"*^) should, of 
course, give the correct answer. In resorting to the saddle-point approximation in real time, 
however, we would have unphysical results with the original factor. It will turn out that 
the above replacement also reproduces the same total tunneling rate as is obtained by the 
instanton method for arbitrary loT^ (see Sec. |III D| ) as well as the same transmission spec- 
trum that is obtained by the transfer-Hamiltonian method [jl6[ in the limit of cuTq (see 

Sec. fvg). 
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The stationary-phase trajectories for X and Y are determined so as to make S{t) ex- 
tremal From conditions 6S{t)/6X{ti) = and 6S{t)/6Y(ti) = 0, we obtain the equations 
of motion for X and Y: 



MXi- 



dV 



-fxX{ti)+i 



J{uj) 



Tx 



MYr 



dX{ti) ' ' Jo n [JO 

Jo TT UJ 

dV 



Ty 



dY{tr 

roo, 



■fiY{ti)+t —J{uj) 



-d 



duj J{uj) 



TT 



JO Jo 



-iul[{tl-TY)-{t2-Tx 



Tl UO 



(32) 



To uniquely determine the solutions to Eqs. (|32D , we must specify initial velocities X{Q) 
and 1^(0). Multiplying the first equation in by X{ti) and integrating it from ti = to 
ti = t', we obtain obtain the law of energy conservation: 



where 



^x'{o) + v{0) = ^x\t') + u^4x{t% 



U,s[Xit')] = V[Xit')] + ^X\t') - I f'dt, r—Jiiu) 

2 Jo Jo IT 

Jo Jo 
r°°dLJ Jiuj) 



(33) 



+d / dtiXiti 



TT UO 



(34) 



is an effective potential which is a local function of time because the stationary-phase trajec- 
tories X and Y are uniquely determined as functions of time. To determine initial velocity 
X(0), we must specify the value of the left-hand side of Eq. (PBD, i.e., the initial energy of 
a tunneling particle. To be consistent with the WKB result, we should identify it with the 



incident energy E of the particle. We thus obtain X(0) = ±iy2[V^(0) — E]/M. To ensure 
that the wave function does not grow but decay inside the barrier, we take the plus sign. 
Initial velocity 1^(0) for trajectory Y can be determined along a similar line, giving 



x(o) = -r(o) 

where the last equation defines Tq. 



l 2[V{0)-E] ^ .d_ 



(35) 



B. Analytic solutions 



The integro-differential equations (p2D can be solved analytically as long as they are 
linear in X and Y . This means that analytic solutions can be obtained if V{X) has the form 
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if X < 0; 

V{X) = ^Uo + U,^ + U2{^y ifO<X<d; (36) 
if X > c/. 

In what follows, we will discuss the case of U2 = to avoid making solutions unnecessarily 
complicated. (The case of f/2 7^ can be treated along a similar line to what is described 
below.) The type of the potential with U2 = appears, for instance, when electrons tunnel 
through a square potential of height Uq which is tilted by voltage bias V, i.e., Ui = —eV. 
Integrating Eqs. (|32D twice with respect to time yields 



X(ti) .ti Ml n 2 f*\ , f°^duj Jiuj) . , 
a Iq 2 MJo Jo TT uj'^ 

I ^ ^V"*i-zcuti-l)a^(t), 



M Jo IT UJ^ 

Y(ti) .ti Ml 2 2 [\ , r°^duj Jioj) . , ,. 

rf!i£ZH(e--.+,.t,-l)6„(t). (37) 

JO TT UJ-^ 



d 



where ui = Ui/Md and 



JO a JO d ' 



b^{t) = -luj dt'e'^''^ dtV-(*'-*-^-+^-)^ + (l-e-*-*)e*^^-. (38) 

JO a JO d 

From Eqs. (|37|) and (^) we find that the solutions have the following general properties: 

x(ti) =r(-ti), 

Tx = -Ty (39) 

These properties can be used to reduce the coupled integro-differential equations to a single 
one, i.e., the first equation of Eqs. (^), which is solved by the Laplace transformation. The 
Laplace transform X[s) of X{t') is given by 

X(,) _ ^.+Mi+^/o-^:y^a.(t) ^ ^^^^ 



S3 



^ + M JO 



M JO TT UJ S^+Ul'^ 



where 



T 



X{t' 



a^it) = luj dt' [e"*"* + e-*"(*+* -2^)J + (1 - e-^"*)e-*"^. (41) 

Inverse Laplace transformation of Eq. (|iO|) gives solution X{t') which includes a^{t) and T 
as undetermined. Here a^{t) is determined by substitution of the obtained solution into 
Eq. (^rp and T is determined from condition X{T) = d. Equations (|39D-(pD therefore give 
the exact solutions for the stationary-phase trajectories for an arbitrary spectral density 
function of the bath. 
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C. Tunneling exponent 



Integrating by parts the kinetic energy terms in Eq. ( |50D and eliminating the resulting 
second time derivatives of X and Y using Eqs. (p2D, we find that 



S{t). 



Mdy,: 



X{Tx)-Y{Ty) +-/ dh 
-I 2 JO 



1 r^^. 



dV 

X^-2V{X) 



dX 



1 pTy 

- / dti 



dV 

Y^-2V{Y) 



. , f^duj J{uo) 
+ldl 



-iujt\ 



X{t, 



10 TC UJ^ 

Relations (|39|) can be used to simphfy this as 

So 2d 2To SoJo ^ d 2mJo tt lu^ 
where 



d 



Jo oj^ Jo a 



5o = 4(f/o - E)To 

is the tunneling exponent in the absence of the interaction. 

Equations (|iO|), PID and ( ^3| ) constitutes the complete solution to the problem. 



(42) 

, (43) 
(44) 



D. Relation to the instanton method 



The effective action S'efr(t) can be related to the effective action obtained by the instanton 
method. The instanton method imposes the boundary conditions X(0) = X{2T) = 0. To 
satisfy them, we replace time variables ti — T {i = 1, 2) by — i(rj — T) for trajectory X(ti) 
and by i{Ti — T) for trajectory Y{ti). The resulting effective action, which we denote as 
S^{t), reads 



5^(t) 



2T 

dr 



M 



+ V{X) + 



r 



""duj 



+d' {nj,{~uj){l 
Jo 2tt 



-iujt 



uj I Tie 
Jo 



2T r-2T 

dTi / (ir2 a(ri - T2)X{ti)X{t2) 
JO 

^^^'^^' + (^^-c^)^ (45) 



i^(ti-T); 



where a(r) = a^{t = —i\T\) and the second term in the curly braces can be obtained by 
replacing u in the preceding term by —u. As shown in Eq. (|26|) , the effective action that 
gives the total tunneling rate is given by setting t = in S^{t = 0). Then the last term in 
Eq. (|45|) , which is responsible for decomposing the total escape rate into transfer energies, 
vanishes and the remaining action takes the form similar to the one obtained by Caldeira 
and Leggett: [|l^ 



S' 



CL 



f3n 
dr 



M 



X^ + V{X) 



+ 



dn / dT2a{n-T2)[X{n)-X{r2)y 

Jo 



(46) 



Different upper bounds between 5*^ and 5""^ result from different boundary conditions. In 
the former case, the boundary condition X(0) = X{2T) is imposed, while in the latter case 
the thermodynamic boundary condition X{0) = X{hf3) is imposed. 

In the following two sections we will use these results to compute the transmission spectra 
for a bath with single frequency and for a bath with Ohmic spectrum. 
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IV. TRANSMISSION SPECTRUM FOR A BATH WITH SINGLE FREQUENCY 



The transmission spectrum for a bath with single frequency has been studied for the case 
with and without the counter term in Refs. fl^ and [1^ using an imaginary-time approach. 
This section has two aims: One is to confirm that the resuhs obtained by the imaginary-time 
approach are reproduced by the real-time approach presented in the preceding sections. 
The other is to derive a new analytic expression of the tunneling exponent which covers from 
sudden tunneling {ujTq ^ 1) to adiabatic one {ujTq ^ 1). These two regimes were treated 
separately in Refs. [p!7| . 



A. Stationary-phase solution 

The spectral density function of a bath with single frequency u is given by 

where 7 is a dimensionless coupling constant. Substituting Eq. p7|) into Eq. 
inversely Laplace-transforming it gives 



(47) 
and 



d 



i ( sinpcutA Ml / 9 . 1 — cosncuti 



+7a(^(t) 1 — cospa;ti -|- iuti — i 



sin pujt I 
P 



(48) 



where p = y/1 + 27. Substituting Eq. ( ^8|) into Eq. (|T|) and solving for a^{t), we find that 



1 2{l-cospujT) + ^{smpujT-pujT)-F{uj){l-e- 



■iuit\ 



ojTq cos pujT +ip sin pujT+2'y — ^ {sin pujT—pujT){l—e 



(49) 



where 



N 1 m , -sinpuT 2 ^ 

r [uj) = l — cospiLiI +t tujl +pujl 



p 



2Ui 



P^ujSo 



1 — ^ ^ cospci;T+ ip sin pooT — ip^ujT 



(50) 



B. Small-7 expansion 

When 7 <C 1, it is sufficient to keep terms up to the first order in 7. Neglecting terms of 
the order of 7^ and those of the order of -yUi/Uo, we find that 



X{t,) 



To 



:i-27)ti + 27 



sin uti 



Ui 



t\ + 'yauj{t) {l — cosuti+iuti — i sinujti). 



(51) 
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Since the last term in Eq. ([51| ) is already multiplied by 7, it is sufficient to expand Eq. 
up to the zeroth orde in 7: 



UjTn 



2(coshcjro - 1) + (1 - e-'^^«)(l - e 



(52) 



To evaluate the tunneling exponent in Eq. (^), we need to find T. Since 7 is small, it is 
reasonable to assume that T/Tq can be expanded in powers of 7: 



T 



-iil + x) 



UiTq 7 



5*0 loTq 



2 (e"'^^"- 1 + cuTo ) - ( 1 - e""^'' ° ) ' ( 1 + cuTo) 



+ ( 1 - e-"^° ) ( 1 - e-'^^°- cuToe-'^^° ) ( 1 - e"^^*)] } . 
Substituting Eqs. into Eq. (|D yields 



where 



5o 



Dfx) 



f/iTo 



5n 



^-iD{uTo) +7iy (cuTo) (1 - e- 



and 



W{x) 



:i-e 



2x 



1-^ ifx>l, 



ifa;<l; 
if X > 1. 



2x 



Substituting Eq. ( pSj ) into Eq. (^) yields the desired transmission spectrum 



r(E,E')=Arexp|-^ 



1 



oj 



n=0 



(53) 



where x is of the order of 7 or Ui/Uq. From condition X{T) = d and Eqs. (|5TD-(p3D, we find 
that 



(54) 



(55) 



(56) 



(57) 



(58) 



Here the first exponent gives the total tunneling rate while the remaining term shows how the 
total tunneling rate is distributed over final states having different energies. Each channel 
n corresponds to a process in which n quanta of the bath are excited by the tunneling 
process. The first exponent on the right-hand side of Eq. (^) comprises three ingredients in 
the square bracket: 1 gives the WKB exponent in the absence of the interaction, —UiTq/Sq 
gives the contribution from the bias, and '~fD{ujTo) describes the effect of the interaction. The 
positive sign of D{ujTo) implies that the total tunneling rate is suppressed by the interaction. 
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Since D{ujTq) becomes unity in the limit of ujTq oo, the effect of the interaction does 
not appear to vanish in this hmit, although the original Hamiltonian (P) appears to require 
that it should vanish because of the presence of the counter term 

/^counter = E 7^X' = J^lduJ. (59) 



Since the counter term gives the maximum possible interaction energy, the ratio of it to the 
bare potential Uq may be regarded as the "physical" coupling constant g which is to be held 
constant. Noting that the counter term scales as (X/ci)^, we are led to define g as 



9 = 



He, 



For the spectral density function in Eq. (E^ we have 



(60) 



« = or (61) 

If we take the limit uTq — >■ oo with g held constant, the effect of the interaction vanishes as 
~ (I/cuTq)^ in the limit of uTq oo, whereas it remains finite if g is increased in proportion 
to uTq as seen above. It is interesting, however, to observe that the single parameter that 
characterizes the tunneling exponent is 7 and not g. In the physics of small-capacitance 
tunnel junctions, g^ gives the ratio of the elementary charging energy to the bare barrier 



height |10]. Figure 1 shows the normalized tunneling exponent S/Sq as a function of uTq for 
several values of g, where 5* = Scs{t = 0). The total tunneling rate decreases with increasing 
coupling constant g or decreasing ujTq. 



C. Tunneling through a fluctuating barrier: classical vs. quantal modulations 

Quantum tunneling through a fluctuating barrier has seen a remarkable resurgence of 
interest since Biittiker and Landauer |]13[ reconsidered the problem of tunneling through a 



barrier that is time-modulated in a prescribed manner. For an opaque barrier, they found 
that the upper and lower sideband intensities, /+ and /_, of transmitted particles having 
incident energy E satisfy the following relation: 

— =tanhc<jTo, (62) 



where u is the frequency of the modulation and Tq is given by 



To = t^^^=- (63) 



For ijjTq <^ 1 a tunneling particle will see an instantaneous potential, while for uoTq ^ 1 it 
will see a time-averaged one. One can therefore expect a crossover between the two regimes 
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around ujTq = 1. In fact, as the frequency of the modulation increases, the ratio (^) rapidly 
increases from to 1 around uTq = 1, which suggests that a time scale for the traversal time 
for tunneling is given by the inverse characteristic frequency at which the crossover occurs. 

The crossover argument provided by Biittiker and Landauer is relevant to many prob- 
lems of physical interest such as those of the dynamic image potential |T8| , p^ ,|3|j20|| , those of 



tunneling through a Josephson junction |21|, and more recently those of tunneling through 
smal-capacitance normal tunnel junctions that could operate at room temperature |TD|,|22[. 
In these problems, however, barrier modulations are caused not by classical sources but by 
dynamical degrees of freedom contributed mostly from their zero-point fluctuations. In the 
problem of dynamic image potential, for example, the barrier modulation is caused by the 
zero-point fluctuations of surface plasmon modes [^. In the problems of tunneling through 
small-capacitance normal and superconducting junctions, barrier modulations are caused 
by quantum fluctuations of the electromagnetic (EM) environment surrounding the junc- 
tion piJTOpi]. We use the results obtained thus far to study the case in which the barrier 



modulation is caused by the quantum-mechanical and thermal fluctuations of the environ- 
mental degrees of freedom. When the coupling of the particle to the bath is not strong, we 
find from Eqs. (^) and (P) that 



sinh^ 

n=— oo L^iiJ-iJ- 2 J 



■(^^/^)S{E' - E -nhu), (64) 



where /„ is the modified Bessel function, and h{Lu) = {'jSo/fiLuTo) sinh^(ti;To/2). For the n-th 
upper and lower sidebands /±„ we find that 



tanhno; Tq . (65) 



In the high-temperature limit where /5 — 0, Eq. (^) reduces to the result obtained for 



classical modulation, esp., a specific result for n = 1 coincides with that obtained in Ref. [13 
In the low-temperature limit where /3 — > cxd, however, Eq. (|65D approaches -1. Physically this 
is because the particle cannot absorb energy from the bath at zero temperature. Thus the 
crossover argument for the traversal time for tunneling does not hold true when the barrier 
modulation is caused mainly by quantum-mechanical zero-point fluctuations. A special case 
of n = 1 is illustrated in Fig. 2 for various values of the ratio of the energy quantum huo of 
the bath to the thermal energy f3~^. 



V. TRANSMISSION SPECTRUM FOR A BATH WITH OHMIC SPECTRUM 

A. Stationary-phase solution 

The spectral density function for a bath with Ohmic spectrum is given by 

J{uj) = Mr]uj. (66) 
Substituting Eq. ( |BUD into Eq. (^) and inversely Laplace-transforming the result yields 
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^ = ^(l-e-''*^) + ^(e-^*^-l + r/tO 



+ —aUt) 

Jo TT 



17] 



uj uj[ri + luj) 



Substituting this into Eq. (|4l| ) and solving for at^(t) gives 

i 



7 ?/ 1 



+i 



'duj'al{t) 



7]^ 
V 



7] + iuo' 



10 TT UJ 

where fu,9i^,K, and are given by 

Auit' _|_ -iujt iui{t' -2T) 



T 



iuj I dt' 
Jo 



with k{t') = 1, e~^^' , e*'^*', and t', respectively. 



(67) 



(68) 



(69) 



B. Small- ?7 To expansion 

When 7]Tq <^ 1, it is sufficient to keep terms up to the ffist order in 7]Tq. Neglecting the 
higher-order terms, we obtain 



X{t 



i) ^ / V 0^ ui r, r°° du a^(t) ,^ . ■ , 

To 2 2 JO TT cu^ 



d 



where 



UoTr 



{I - e-'^^){l - e 



-iujT\ ^—iujT^—iujt'^ 



(70) 



(71) 



The behavior of Eq. ([70|) at large t determines the transmission spectrum near zero bias, 
while its value at t = determines the total tunneling rate. From condition X{T) = d we 
find that 



T 



1 + ^1^0 + 



(72) 



So ' 2n 

We ffist consider the total tunneling rate which is determined by S{t = 0). We find that 



S{t = 0) 



l-^ + ^ln2 



(73) 



The positive sign of the last term in the square bracket indicates that the total tunneling 
rate is suppressed by the interaction with the Ohmic bath p!2 |. 

We next consider the transmission spectrum near zero bias E E' which is dominated 
by the value of S{t) at large t: 
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So 



(74) 



Substituting this into Eq. p9| ) we find that 
T{E,E')=N 



exp 


_5a 


(1 - ^ + l^^^o) 


1 


[i? — E'\Tq 




|E - E'\ 


h 



(75) 



C. Zero-bias anomaly 

Equation (|75D imphes that the differential conductance obeys a power law near zero bias. 
To be specific, we consider the case of electron tunneling between electrodes having the same 
chemical potential, and assumes for simplicity that the density of states, Dl, -Dr, of the 
left and right electrodes are constant. The tunneling current liV) at bias voltage V is then 
given by 



I{V)=eD^D^JdE JdE'T{E, E-E') {f{E) [1- f{E~E' + eV)] 
-f{E) [l-f{E~E'-eV)]}. 



At zero temperature, this reduces to 



I[V) = eDj^D^ r^dET{E, E-E'){eV - E') 
Jo 

Substituting Eq. (|75|) into this yields 



(76) 



(77) 



IiV)=e^DLDKN 



exp 



-f (1 



UiTq 



So +4;^^^0 



2 + r (^r^To) 



V 



eVTo 



So 



vTo 



(78) 



Here is determined by the requirement that the differential tunneling resistance in the 
absence of the interaction at zero voltage is given by Rt'- 



(79) 



We thus find that 



exp 


'^{eVTo-lvToSo)' 


V 


\eVTol 






r (2 + ^,vTo) 




[ n \ 





(80) 



This indicates that the differential conductance at low voltages V obey a power law 
(el^To//^)''^o^«/2.^ 
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Fig. 1. Normalized tunneling exponent S/Sq as a function of uoTq for three values of g, 
where S = Scs{t = 0), and 5*0 and Tq are the tunneling exponent and the barrier traversal 
time in the absence of interaction. 

Fig. 2. Ratio of the first sideband intensities I± = I±i plotted as a function of uTq for 
various values of j3hu, where is the thermal energy and hu; is the energy quantum of 
the single-frequency bath. 
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